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Introduction 


Oxygen  toxicity  is  a  consequence  of  breathing  oxygen _ partial 
pressures  greater  than  in  atinospheric  air.  While  oxygen  itself 
is  not  toxic,  the  evidence  indicates  that  toxic  derivatives  of 
oxygen  are  a  by-product  of  cellular  respiration.  The  production 
of  toxic  oxygen  species  increases  as  the  tissue  oxygen  tension 
rises  (Yusa  et  al.  1987),  and  toxicity  probably  occurs  due  to 
tissue  damage  or  interference  with  normal  function.  The ^ toxic 
oxygen  species  produced  in  normoxic  conditions  are  deactivated, 
but  under  hyperoxic  conditions  they _ accumulate.  The  balance  be- 
production  and  removal  of  toxic  species  is  well— suited  to 
mathematical  modelling  and  statistical  analysis . 

Gerschman  (1964)  was  the  first  to  suggest  that  oxygen 
toxicity  might  be  caused  by  free  radicals  which  are  highly  reac¬ 
tive  and  potentially  damaging  to  cellular  components.  Although 
the  valididty  of  this  hypothesis  remains  to  be  prcpven  absolutely, 
evidence  has  accumulated  associating  oxygen  toxicity  with  in 
creased  generation  rates  of  reactive  oxygen  species  such  as  su¬ 
peroxide  (O2-),  hydrogen  peroxide  (H2O2),  singlet  oxygen,  and 
hydroxyl  radicals  (Freeman  and  Crapo  1981). 

The  principal  reactions  leading  to  these  oxygen  species  are 
(Fridovich  1979;  Freeman  and  Crapo  1981). 

(O2)  +  electron ^  (02“ ) 

SOD 

2  (O2-)  +  2  (H+)  (H2O2)  +  (O2) 

CAT 

2  ( H2O2 )  (^2^  ^  (H2O) 


where  the  second  and  third  reactions  are  catalyzed  by  superoxide 
dismutase  (SOD)  and  catalase  (CAT)._  If  these  reactions  do  not 
pj-oceed  to  completion,  hydroxyl  radicals  can  also  be  formed. 

The  formation  of  a  small  number  of  reactive  oxygen  species 
is  a  normal  by-product  of  cellular  respiration  (Fridovich  1979). 
Most  of  these  are  deactivated  by  SOD  and  CAT  or  may  be  scavenged 
by  vitamins  A,  C,  and  E.  Unscavenged  oxygen  species  can  damage 
cell  membranes  and  intracellular  enzymes  by  reacting  with  tissue 
proteins  and  lipids .  This  damage  is  repaired  by  cellular  defense 
and  biosynthetic  mechanisms  so  that  acute  oxygen  toxicity  does 
not  occur  under  normoxic  conditions.  If  hyperoxia  increases  the 
production  of  reactive  species  beyond  the  deactivating 
capabilities  of  SOD  and  CAT  and  beyond  the  capacity  of  the 
reparative  mechanisms ,  then  cellular  damage  accumulates  and 
eventually  leads  to  symptoms  of  oxygen  toxicity  (Clark  1982; 
Fisher  et  al  1979). 
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PULMONARY  OXYGEN  TOXICITY 
Clark-Lambertsen  Free  Radical  Model 

The  UPTD  concept  introduced  by  Clark  and  Lambertsen  (1970) 
to  track  the  development  of  pulmonary  oxygen  toxicity  may  be  in¬ 
terpreted  as  a  balance  between  production  and  removal  of  a  toxic 
substance.  The  UPTD  concept  was  based  upon  a  graphical  analysis 
of  decreases  in  vital  capacity  which  occur  during  the  inspiration 
of  hyperoxic  gas  mixtures.  These  decreases  were  described  by 
.rectangular  hyperbolas  having  the  form 

(PO2  -  0.5  ATM)t  =  k 

where  0.5  ATM  is  the  estimated  threshold  PO2  above  which  pul¬ 
monary  toxicity  begins  to  occur,  and  the  constant  k  is  an  index 
of  the  decrease  in  vital  capacity. 


This  relationship  may  be  derived  theoretically  from  a  mass 
balance  for  a  toxic  substance  (F).  Consider  a  tissue  in  which  F 
is  produced  at  a  rate  proportional  to  the  local  PO2  and  removed 
at  a  constant  rate.  The  mass  balance  for  F  gives 


dF/dt  =  C2  PO2 


=  2 


Integrating  F  with  respect  to  time  (t), 

F(t)  =  Fq  +  (Cl  PO2  -  C2)t  (1) 

where  Fq  is  the  initial  concentration  of  F  and  and  C2  are  con¬ 
stants.  This  can  be  re-arranged  into  the  form  given  by  Clark  and 
Lambertsen  (1970). 


(PO2  -  C2/C1)  =  (F  -  Fq)/ci 


(2) 


Neglecting  the  exponent  of  t,  which  Bardin  and  Lambertsen ^ 
(1970)  later  raised  to  1.2,  the  Pulmonary  Toxicity  Dose  (PTD)  is 
defined  as 

PTD  =  C3  (PO2  -  C2/ci)t  (3) 

where  C3  =  2  and  ~  ATM. 

Rearranging  these  terms  gives  the  familiar  expression 

PTD  =  t(0.5/(PO2-0.5) )~^ 


and  the  percent  change  in  vital  capacity  (VC)  is  approximately 

VC  =  -  C4  (PTD  -  C5) 
where  C4  =  0.01  and  C5  =  415. 


(4) 
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This  formulation  of  the  Clark-Lambertsen  concept  ^nay 
rrivpn  the  following  interpretation.  Under  normoxic  condition  , 
a::?  Ses  n^  ha5I  pulmoLry  oxygen  toxicity  and  his  vital  capa- 
unaffected.  According  to  equation  (1),  all  free  radicals 

are  removed  and  F  =  0  since 

Cl  PO2  <  C2 

Thus,  the  initial  free  radical  concentration  is  zero  or 


By 

constants 


Fo  =  0 

combining  equations  (2)  -  (4)  and  eliminating  redundant 

Vc  =  -  C4  (F  -  C5 )  ( 5  ^ 

Eauations  (1)  and  (5)  define  a  model  which  can  follow  changes  in 
vital  capacity  during  any  oxygen  partial _ pressure  profile.  Whe 
experimental  measurements  of  vital  capacity  are  availabe,  more 
over  the  error  between  measured  and  estimated  -s  VC  change  can  b 
determined  and  the  best  values  for  the  model  parameters  can  be 
found  by  minimizing  this  error. 

Vital  Capacity  Data 

The  Clark-Lambertsen  model  was  applied  to  the  vital  capacity 
data  compiled  by  Harabin  et  al  (1986).  All  control  measurements 
before  the  start  of  hyperoxic  exposure  were  included  in  the  data 
set  diving  794  vital  capacity  measurements  rather  than  the  -40 
Itudie^S?  Harabin  et  al  (1986).  Data  points  for  each  subDect 
were  converted  from  liters  to  percent  vital  capacity  change  by 
dividing  by  the  first  control  measurement  for  each  subject. 

Vital  Capacity  at  Low  PIO:) 

The  use  of  vital  capacity  as  an  index  of  pulmonary  oxygen 
toxicity  f?Luently  criticized  because  of  variability  in 
rSeated  mlasSSments!  This  is  reflected  in  Fig.  1  which  is  a 
freauency  diagram  of  312  vital  capacity  measurements  made  at 
PIOo's  of  0.21  to  0.5  ATM  before  exposure  to  higher  PIO2 
percent  change  in  vital  capacity  was  tabulated  over  increments  o 
2°-.  The  mean  PIOo  and  %  VC  change  of  this  distribution  are  0*296 
ATM  and  -1  11%.  Thus,  use  of  the  first  vital  capacity  measur 
ment “o  djilne  the  normal  vital  capacity  for  each  subject  un- 
derestimated  the  mean  by  1%. 

Ninety-five  percent  of  this  distribution  (the  ^ 

Panap")  falls  between  %  VC  changes  of  -11  and  +8^.  This  suggesps 
?ha?  a  vi?al  capacity  measurement  must  differ  from  the  i^ean  by 
SSoximitely  ±9%  fo?  there  to  be  a  95%  chance  that  the  measure- 
ment  is  valid  and  not  a  statistical  aberation. 
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Figure  2  subdivides  the  measurements  from  Fig.  1  into  3 
smaller  ranges  as  shown  in  Table  1. 

Table  1.  Percent  change  in  vital  capacity  at  PI02's  of  not  more 
than  0 . 5  ATM 


No.  of 

Measurements 

PO2  Range 
(ATM) 

Mean  PO2 
(ATM) 

Mean  %  VC 
Change 

95%  VC 
Range 

74 

0.21-0.23 

0.214 

-0.924  ■ 

-13%  to  +7% 

130 

0.28-0.34 

0.298 

-1.117 

-  9%  to  +5% 

48 

0.47-0.5 

0.482 

-2.263 

-13%  to  +4% 

312 

0.21-0.5 

0.296 

-1.110 

-11%  to  +8% 

The  number  of  measurements  in  the  first  three  rows  does  not 
To  then^er  ?n  the  fourth  row  to  avoid  the  influence 
exist)  of  0.5  ATM  PIO2 '  s  on  subsequent  VC 
PI02's.  The  95%  VC  Range  varies  from  14^  to  200  for  the 

PIO2  sub-ranges. 

While  there  are  no  striking  differences  the  3  ^istrib- 
ii^inns  of  Fia  2  Fig.  3  shows  a  plot  of  the  mean  WC  change 
aaa?Slt  the  rakn  PIO,  for  each  distribution  and  suggests  a  decre- 

^St  in  capacity  With  increasing  PlOg  at  ualues^of^up  to 

n  ^  ATM  Tt  is  coininorilv  3.ssTjjn0ci  tlri3.ti  3.  *  -i-  x. 

ttLrated  with  nSShangS  in  vital  capacity,  and  this  may  be  true 
toierat  Harabin  et  al  (1986),  however,  demonstrated 

^?anificant  variation  in  individual  response  to  elevated  oxygen 
?a??iifSessure,  and  susceptible  individuals  may  be  responsible 
for  the  decrements  in  Fig.  3. 


Vital  Capacity  Data  Analysis 

The  vital  capacity  data  was  ^"alyzed  by  finding  the  best  fit 

^aii  “?nS"VhJ^  fl/was 

accomplished 'by  finding  the  values  of  the  model  parameters 
accomplisnea  py  r  J  -^^3  d)  and  (5))  which  minimized  the 

errors  (SSE)  between  the  predictions  of  the  _ 

Sel  and  the  data.  The  minimization 

=trsss^irs^r^rrLgre=ft^rsttrfofLcruS2tii^i^id  para- 

meter. 

While  manv  free  radical  and  vital  capacity  models  were 
cal  model  described  earlier  and  designated  FR(1),  another  tre 
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radical  model  (FR{2))  assumed  that  free  radical  elimination  had  a 
t6riTi  proportiondl  to  its  own  concontmtion  or 

dF/dt  =  Cl  PO2  ”  {C2  +  C3  F) 

This  equation  has  the  solution, 

F(t)  =  (Fq  -  C4)  e-'^^^t  +  C4 

where 


C4  (C]_  PO2  -  C2)/C3. 

The  vital  capacity  model  of  eg^s-tion  (5)  was  simplified  to 


VC  =  -F 


in  which  each  rectangular  hyperbola  defined  by  equation  ( 1 )  cor 
responds  to  a  different  decrement  in  vital  capacity.  This  is 
designated  as  vital  capacity  model  VC(1). 


Model  VC(1)  predicts  that  the  vital  capacity  can  decrease  to 
less  than  -100%.  -  As  this  is  unrealistic,  model  VC(2)  was 
defined  in  which  the  vital  capacity  has  a  minimum  limiting  value 
of  -100%.  Thus, 

VC  =  -100  ( 1  -  b3_“^)  . 


The  free  radical  and  vital  capacity  models  are  summarized  in 
Table  2. 


Table  2.  Principal  free  radical  and  vital  capacity  models. 


- [ 

N 

FR(N) -Free  Radical  Model 

VC (N) -Vital  Capacity  Model 

1 

dF/dt  =  Cl  PO2  -  c  2, 

< 

0 

11 

1 

F(t)  =  Fq  +  (Cl  PO2  -  C2)  t 

2 

dF/dt  =  Cl  PO2  -  (C2  +  C3  F) 

F(t)  =(Fo-C4)  e  +  C4 

VC  =  -100  (1-bi^^ 

C4  =  (Cl  PO2  -  C2)/C3 

For  an  analysis  of  vital  capacity  data,  the  PO2  in  the 
models  of  Table  2  is  the  alveolar  value  in  the  lungs.  This 
was  defined  as 


value 
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PA02  =  PI02  -  PH2O  -  PAC02 
where  PH2O  =  47  torr  and  PACO2  =  40  torr 

PAO2  =  PIO2  -  0,053  ATM. 

Model  FR(1)-VC(1)  is  identical  to  models  used  by  Harabin  et 
al  (1986).  In  their  analysis,  vital  capacity  measurements  were 
treated  as  discrete  points.  In  the  analysis  reported  here,  each 
vital  capacity  measurement  depends  upon  the  previous  measure¬ 
ment.  Thus,  decreases  in  vital  capacity  develop  and  resolve  con¬ 
tinuously  for  individual  subjects.  The  Harabin  report  analyzed 
440  measurements  while  this  report  analyzed  794.  These  included 
all  measurements  made  at  POl2's  of  0.5  ATM  or  less. 

Table  3.  Results  of  dataset  analysis  with  various  models.  Add 
0.053  ATM  to  C2/c;j_  to  arrive  at  the  PIO2  which  can  be 
tolerated  indefinitely.  More  exact  parameter  values 
are  given  in  Appendix  A.  The  numbers  in  parentheses 
are  from  Harabin  et  al  (1986). 


Free  Radical 
Model  & 
Constants 

— 

Vital  Capacity 
Model  & 
Constants 

— 

SSE 

TT5 

TT6 

1 

FR(1) 

Cl  =  0.011 

C2/ci=0.5  atm 

VC  ( 1 ) 

145,179 

-2.20 

(-2.3) 

2 

FR(1) 

C2^=0 . 009 
C2/ci=0.38  atm 

VC(1) 

156,427 

-1.95 

(-2.0) 

-3.79 

(-3.0) 

3 

FR(1) 

0^=0.0024 
C2_/c2=0.28  atm 

VC(1) 

55,948 

-0.56 

-1.09 

4 

FR(1) 

C2=0.021 
C2/ci=0.41  atm 

VC(2) 

bi  =  1.0055 

31,538 

-2.43 

-4.66 

5 

FR(2) 

C]_=0.0092 
C2/ci=0.17  atm 
C3=0. 00031 

VC(1) 

37,214 

-2.21 

-4.23 

6 

FR(  2  ) 

C2=0.021 
C2/c2=0.41  atm 
C3=-0. 28x10“/ 

VC(2) 

bi=1.0055 

31,518 

-2.86 

-5.05 
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Table  3  presents  the  model  parameter  values,  Sura  of  Squared 
Errors  (SSE),  and  predicted  vital  capacity  changes  for  Oxygen 
Treatment  Tables  5  and  6  (TT5  and  TT6).  Rows  1  and  2  of  Table  3 
use  the  "minimized"  and  "simplified"  parameters  from  the  Harabin 
analysis,  and  Row  3  presents  the  results  for  the  same  model  op¬ 
timized  as  model  FR(1)-VC(1).  There  are  large  differences  in  SSE 
between  the  Harabin  models  (Rows  1  and  2)  and  the  optimized  model 
(Row  3)  which  is  not  unexpected  given  the  different  datasets  and 
applications  of  the  model.  This  points  out  the  importance  of  the 
dataset  in  determining  the  behavior  of  a  given  model.  Indeed,  it 
might  be  said  that  a  dataset  is  represented  by  the  parameter 
values . 

The  %  VC  changes  for  TT5  and  TT6  given  in  parentheses  in 
Rows  1  and  2  are  the  values  computed  by  Harabin  et  al  (1986).  In 
Row  3,  the  value  of  the  oxygen  asymptote  (02/0^  +  0.053)  is  0.333 
ATM,  below  the  commonly  used  value  of  0.5  ATM.  The  %  VC  changes 
for  TT5  and  TT6  are  only  half  of  those  predicted  in  Rows  1  and  2. 

Row  4  is  the  optimized  model  FR(1)-VC(2).  Placing  a  limit 
on  the  vital  capacity  decrement  substantially  improves  the  SSE 
and  increases  the  oxygen  asymptote  to  0.463  ATM.  Row  5  is  model 
FR(2)-VC(1).  The  SSE  is  improved  over  FR(1)-VC(1),  but  the 
oxygen  asymptote  is  now  0.223  ATM  which  seems  unreasonably  low. 
Comparing  FR(1)-VC(2)  and  FR(2)-VC(2)  in  Rows  4  and  6,  there  is 
little  difference  in  SSE,  and  the  value  of  C3  is  small.  Thus, 
there  appears  to  be  little  justification  for  using  a  model  more 
complicated  than  the  modified  version  of  Clark-Lambertsen  desig¬ 
nated  FR(1)-VC(2) . 

Among  other  models  tested  were  sigmoidal  vital  capacity  re¬ 
lationships,  an  F*P02  term  in  the  free  radical  model,  and  a  vital 
capacity  model  incorporating  the  time  integral  of  F.  None  of 
these  models  performed  as  well  as  FR(1)-VC(2). 

Recoverv  of  Vital  Capacitv  Decrements 


One  of  the  more  interesting  aspects  of  oxygen  toxicity  mod¬ 
elling  is  the  recovery  of  vital  capacity  decrements  upon  return 
to  near  normoxic  breathing  gas.  Figure  4  is  the  oxygen  partial 
pressure  profile  described  by  Eckenhoff  and  Vann  (1985)  during 
decompression  from  a  saturation  dive  which  tested  the  development 
and  recovery  from  pulmonary  oxygen  toxicity  in  12  divers  as  indi¬ 
cated  by  changes  in  vital  capacity.  The  vital  capacity  measure¬ 
ments  were  reported  by  Harabin  et  al  (1986).  Figure  5  is  taken 
from  Harabin 's  report  and  is  a  graphical  representation  of  the  % 
VC  changes  for  all  12  divers. 

Figure  6  is  the  prediction,  using  model  FR(2)-VC(2)  from  Row 
6  of  Table  3,  of  the  %  vital  capacity  changes  during  the  oxygen 
partial  pressure  profile  of  Fig.  4.  The  mean  vital  capacity 
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changes  for  the  12  divers  from  Fig.  5  are  also  shown  in  Fig.  6. 
The  predicted  and  mean  profiles  agree  quite  well  during  the 
decrease  in  vital  capacity  while  the  PIO2  is  1.05  ATM.  There  is 
essentially  no  change  in  predicted  VC  when  the  PIO2  is  0.5  ATM, 
but  the  experimental  mean  may  undergo  a  decrease  indicating  some 
recovery.  The  fluctuations  during  this  period  may  reflect  the 
individual  variabiity  which  Harabin  et  al  (1986)  noted  in  therr 
analysis.  During  the  subsequent  recovery  phase,  the  predicted _ 
and  mean  values  are  in  good  agreement  except  for  the  final  period 
in  which  predicted  recovery  .is  faster  than  the  experimental  mean. 

It  is  also  of  interest  to  note  that  PIO2  =  0.59  ATM  in  the  5 
min  air  breaks  during  TT5  and  TT6.  Thus,  these  air  breaks  may 
provide  no  recovery  from  pulmonary  toxicity. 

CNS  OXYGEN  TOXICITY 


Cerebral  Oxygen  Tension 

The  free  radical  models  developed  above  also  can  be  applied 
to  an  analysis  of  CNS  toxicity  data.  In  the  case  of  CNS 
toxicity,  the  relevant  tissue  PO2  is  substantially  lower  in  the 
brain  than  in  the  lungs  for  pulmonary  toxicity.  To  find  the 
lower  PO2  in  the  brain,  the  following  procedure  was  adopted: 

1.  Determine  the  PAO2  as  described  earlier: 

PAO2  =  PIO2  -  0.053  ATM 

2.  Use  the  Hill  equation  to  relate  the  PAO2  to  the  oxygen  con¬ 
tent  of  arterial  blood  (CAO2).  The  physically  dissolved  oxygen 
must  be  included.  Thus 


Ca02  =  a02*PAO2  +  C100/(l  +  (a^  PaOl ^  ^ 

where  a  O2  =  0.023/ATM  is  the  physical  solubility  of  oxygen  in 
blood , 

ClOO  =  0.197  ml  02/ml  blood  is  the  oxygen  carrying  capacity 
of  saturated  hemoglobin,  and 

aQ_  =  24  and  a2  =  2  are  constants  to  fit  the  Hill  equation  to 
the  Hb02  dissociation  curve. 

3.  Use  the  Fick  relationship  to  find  the  venous  oxygen  content 
(CVO2)  after  assuming  values  for  cerebral  blood  flow  (Q)  and 
cerebral  oxygen  consumption  (m) .  Thus, 

m 

CVO2  =  Ca02  -  Q 
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where  m  and  Q  are  assumed  to  be 

m  =  3  ml  02/min/lOO  g  tissue 
Q  =  55  ml/min/100  g  tissue. 

4.  Use  the  relationship  in  Step  2  to  find  PVO2  from  knowledge  of 
CvOt.  This  relationship  cannot  be  solved  explicitly  for  ^^^2 
must  be  solved  numerically.  Appendix  B  tabulates  this  solution 
which  is  available  in  the  computer  program  as  a  look-up  table. 

The  depth  in  Appendix  B  is  the  water  depth  at  which  the  diver  is 
assumed  to  breathe  100%  oxygen.  The  PIO2,  PAO2,  and  PVO2  are^ 
also  given.  It  is  assumed  that  the  PVO2  is  the  relevant  para 
meter  for  the  development  of  CNS  oxygen  toxicity. 

The  numerical  values  assumed  for  the  various  constants  in 
the  procedure  described  above  were  allowed  to  vary  as 
dertermined  parameters  when  the  models  were  fit  to  the  data. 
Surprisingly,  they  had  little  influence  on  the  final  result. 

CNS  Oxygen  Toxicity  Data 

The  human  CNS  oxygen  toxicity  data  used  in  the  analysis  was 
from  U.S.  Navy  studies  published  by  Lanphier  (1954),  Butler  an 
Thalmann  (1984;  1986),  Piantadosi  et  al  (1979),  Schwartz  (1984), 
and  Butler  (1986).  Consistent  with  these  reports  and  current 
NEDU  practice,  symptoms  of  CNS  oxygen  toxicity  were  definedas 
Probable,  Definite,  or  Convulsions.  The  data  was  analyzed  in  3 
symptom  categories : 

1.  Convulsions  only  (11  occurrences) 

2.  Convulsions  and  Definite  Symptoms  (44  occurrences) 

3.  Convulsions,  Definite,  and  Probable  Symptoms  (81  occurrences) 
The  total  number  of  man-exposures  was  773. 

Models  and  Analysis 

Symptoms  were  treated  as  a  binary  variable:  Owhen  absent 
and  1  when  present.  This  allowed  the  method  of  maxira^  _ 

likelihood  to  be  used  as  described  by  Vann  (1987)  with  a  variable 
parameter  step  size  introduced  to  accelerate  convergence. 

The  models  FR(1)  and  FR(2)  were  used  to  describe  the 
kinetics  of  the  toxic  substance.  Sigmoidal  probability  models 
were  used  to  relate  the  free  radical  concentration  (F)  to  the 
probability  (P)  of  CNS  toxicity • symptoms .  The  most  successful 
probability  model,  however,  was  similar  to  VC(2), 

P(F)  =  1  -  bi"^. 

Indeed,  while  the  maximum  log  likelihoods  were  somewhat  less  for 
FR(2)  and  the  sigmoidial  probability  models,  some  of  the  op- 
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timized  parameters  were  found  to  have  values  of  zero  whioh 
rendered  the  models  of  no  predictive  use. 

The  most  successful  model  was,  again,  the  simplest:  the 
modified  Clark-Larabertsen  FR(1)  and  the  probability  model  given 
above.  Appendix  C  presents  parameter  values  for  analysis  of  the 
3  symptom  categories. 

Using  the  oxygen  asymptote  defined  by  c-^/c2  and  Appendix  B, 
the  following  predictions  are  made  for  the  depth  on  100%  oxygen 
at  which  a  diver  has  zero  risk  of  developing  CNS  toxicity: 


Svmptom  Category 

Max  Safe 

-02-J 

1. 

Convulsions  only 

19 

fsw 

2. 

Conv.  and  Def.  Sx. 

20 

fsw 

3. 

Conv. ,  Def.,  and  Prob.  Sx 

16 

fsw 

Appendix  C  also  gives  the  probability  of  symptoms  for  the  old  and 
new  U.S.  Navy  oxygen  exposure  limits  and  the  oxygen  exposure 
times  at  5  fsw  depth  increments  for  symptom  probabilities  of  0.01 
to  0-08  in  increments  of  0.01.  For  the  old  and  new  exposure 
limits  at  25  fsw,  for  example,  the  following  probabilities  were 
found: 


Sx  Prob  at  25 

fsw  for 

Svmptom  Category 

75  min 

240  min 

1. 

Conv.  only 

0.008 

0.025 

2. 

Conv.  Sc  Def. 

0.031 

0.095 

3. 

Conv.,  Def.  Sc  Prob 

0.061 

0.183 

An  example  of  the 

development  and  resolution  of 

the  prob- 

ability  of  CNS  toxicity  symptoms  is  shown  in  Figs.  7  and  8.  Fig¬ 
ure  7  shows  an  oxygen  dive  profile  with  120  min  at  20  fsw,  15  min 
at  40  fsw,  and  105  min  at  20  fsw.  Figure  8  shows  the  percent 
risk  of  convulsions  or  definite  symptoms  (Category  2).  This  risk 
increases  slowly  at  20  fsw  to  about  0.2%,  rapidly  at  40  fsw  to 
about  4.1%,  and  slowly  at  20  fsw  to  about  4.3%.  Upon  surfacing, 
the  risk  drops  to  0%  in  10  minutes. 

CONCLUSIONS 

It  is  evident  that  both  pulmonary  and  CNS  oxygen  toxicity 
can  be  modelled  with  simple  assumptions  and  these  models  fit  to 
existing  bodies  of  data.  The  human  data  currently  available  is 
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limited  but  allows  preliminary  estimates  to  be  made.  These 
should  be  considered  relative  estimates  which  can  be  used  for 
planning  future  experiments  and  perhaps  providing  rough  opera¬ 
tional  guidelines.  With  additional  data  and  better  models,  im¬ 
provements  will  probably  follow.  The  methods  also  should  be 
readily  applicable  to  the  analysis  of  animal  experiments. 
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OXCCtDlU  >0 


APPENDIX  A  -  PARAMETER  VALUES 
FR( 1 )-VC(1 ) 

Cl  =  0. 24254555E-02 
C2/C1  =  0. 27838695EH-00 
SSE  =  55948.39453125 

FR( 1  )-VC(2) 

Cl  =  0. 2 1 00484  IE-0  1 
C2/C1  =  0.406'59472E  +  00 
B1  =  0. 1 0055144E+01 
SSE  =  31538.12695312 

FR(2)-VC(1  ) 

Cl  =  0. 92303315E-02 
C2/C1  =  0. 17013633E+00 
C3  =  0. 30700679E-03 
SSE  =  37213.92187500 

FR(2)-VC( 2) 

Cl  =  0. 21297863E-01 
C2/C1  =  0.41017273E+00 
C3  =-0.28009774E-07 
B1  =  0. 10055192E+01 
SSE  =  31518.07421875 


FOR  PULMONARY  TOXiCITY  MODELS 


APPENDIX  B  -  CEREBRAL  P02  AS  A  FUNCTION  OF  OXYGEN  DEPTH 


DEPTH 

PI02 

PA02 

PV02 

DEPTH 

PI0  2 

PA02 

PV02 

0.  000 

0.  210 

0.  157 

0.  059 

31.  000 

1.  939 

1.  886 

0.  147 

1.  000 

1.  030 

0.  977 

0.  091 

32.  000 

1.  970 

1.  917 

0.,  151 

2.  000 

1.  061 

1 .'  008 

0.  092 

33.  000 

2.  000 

1.  947 

0.  155 

3.  000 

1.  091 

1.  033 

0.  093 

34.  000 

2.  OQO 

1.  977 

0.  159 

4.  000 

1.  121 

1.  068 

0.  095 

35.  000 

2.  061 

2.  008 

0.  163 

5.  000 

1.  152 

1.  099 

0.  096 

36.  000 

2.  091 

2.  038 

0.  167 

6.  000 

1.  132 

1.  129 

0.  097 

37.  000 

2.  121 

2.  068 

0.  172 

7.  000 

1.  212 

1.  159 

0.  099 

38.  000 

2.  152 

2.  099 

0.  178 

a.  000 

1.  242 

1.  189 

0.  100 

39.  000 

2.  182 

2.  129 

0.  183 

9.  000 

1.  273 

1.  220 

0.  101 

40.  000 

2.  212 

2.  159 

0.  188 

10.  000 

1.  303 

1.  250 

0.  103 

41.  000 

2.  242 

2.  189 

0.  195 

11.  000 

1.  333 

1.  280 

0.  104 

42.  000 

2.  273 

2.  220 

0.  201 

12.  000 

1.  364 

1.  311 

0.  105 

43.  000 

2.  3D3 

2.  250 

0.  208 

13.  000 

1.  394 

1.  341 

0.  107 

44.  000 

2.  333 

2.  280 

0.  216 

14.  000 

1.  424 

1.  371 

0.  108 

45.  000 

2.  364 

2.  311 

0.  224 

15.  000 

1.  455 

1.  402 

0.  Ill 

46.  000 

2.  3?4 

2.  341 

0.  233 

16.  000 

1.  485 

1.  432 

0.  112 

47.  000 

2.  424 

2.  371 

0.  243 

17.  000 

1.  515 

1.  462 

0.  113 

48.  000 

2  .  455 

2.  402 

0.  254 

18.  000 

1.  545 

1.  492 

0.  116 

49.  000 

2.  485 

2.  432 

0.  266 

19.  000 

1.  576 

1.  523 

0.  117 

50.  000 

2.  515 

2.  462 

0.  278 

20.  000 

1.  <606 

1.  553 

0.  120 

51.  000 

2.  545 

2.  492 

0.  292 

21.  000 

1.  <636 

1.  583 

0.  121 

52.  000 

2.  576 

2.  523 

0.  307 

22.  000 

1.  667 

1.  614 

0.  124 

53.  000 

2.  606 

2.  553 

0.  322 

23.  000 

1.  697 

1.  644 

0.  125 

54.  000 

2.  636 

2.  583 

0.  338 

24.  000 

1.  727 

1.  674 

0.  128 

55.  000 

2.  667 

2.  614 

0.  357 

25.  000 

1.  753 

1.  705 

0.  130 

56.  000 

2.  697 

2.  644 

0.  375 

26.  000 

1.  733 

1.  735 

0.  133 

57.  000 

2.  727 

2.  674 

0.  396 

27.  000 

1.  313 

1.  765 

0.  136 

58.  000 

2.  753 

2.  705 

0.  417 

28.  000 

1.  343 

1.  795 

0.  138 

59.  000 

2.  788 

2.  735 

0.  438 

29.  000 

1.  879 

1.  826 

0.  141 

60.  000 

2.  818 

2.  765 

0.  462 

30.  000 

1.  909 

1.  356 

0.  145 

appendix  C  -  PARAMETER  VALUES  AND  RISK  PREDICTION  FOR  CNS  TOXICITY  MODELS 

SYMPTOM  CATEGORY  ONE 

2.  0138952732  = 

O. 0000000000 
O.  0117034214  =  Cl 
O.  1 175345778  =  C2/C1 

-  MAXIMUM  LOG  LIKELIHOOD 

NUMBER  OF  SYMPTOMS  =  11. 

SYMPTOM  CATEGORY  =  1 

NUMBER  OF  EXPOSURES  =  773. 

FRACTION  WITH  SYMPTOMS  =  0.014 

MODEL  INDEPENDENT  LIKELIHOOD  =  -9.  8026247 


TIME  AT  INDICATED  DEPTH  FRACTIONAL  RISK  BELOW  LINE 


DEPTH 

0. 01000  C 

1.  02000 

0.  0300  0 

0.  04000 

0.  05000 

0.  06000 

0.  07000 

0.  08000 

10.  00 

-82.  3 

-165.  5 

-249. 

5 

-334.  3 

-420.  1 

-506.  7 

-594.  3 

-682.  9 

15.  00 

-175.  0 

-351.  8 

-530. 

5 

-710.  9 

-893.  3 

-1077.  6 

-1263.  8 

-1452.  1 

20.  00 

557.  0 

1 119.  7 

1683. 

1 

2262.  5 

2842.  8 

3429.  3 

4022.  0 

4621.  2 

25.  00 

96.  4 

193.  7 

292. 

1 

391.  4 

491.  3 

593.  3 

695.  9 

799.  5 

30.  00 

45.  1 

90.  6 

136. 

7 

183.  2 

230.  1 

277.  6 

325.  6 

374.  1 

35.  00 

26.  9 

54.  0 

81. 

5 

109.  2 

137.  2 

165.  5 

194.  1 

223.  1 

40.  00 

17.  4 

34.  9 

52. 

6 

70.  5 

88.  6 

106.  9 

125.  4 

144.  1 

45.  00 

11.6 

23.  2 

35. 

0 

46.  9 

59.  0 

71.  1 

83.  4 

95.  9 

50.  00 

7.  7 

15.  4 

23. 

2 

31.  1 

39.  1 

47.  2 

55.  3 

63.  6 

55.  00 

5.  1 

10.  3 

15. 

6 

20.  8 

26.  2 

31.  6 

37.  1 

42.  6 

60.  00 

3.  6 

7.  2 

10. 

8 

14.  5 

18.  2 

21.  9 

25.  7 

29.  6 

DOSE^^’ 

0.  01436 

0.  02B86 

.  0.  043  51 

0.  05831 

0.  07327 

0.  08838 

0.  10366 

0.  11910 

U.  S.  NAVY 

OLD  LIMITS 

NEW  LI MITS 

DEPTH 

TIME  RISK 

TIME 

RISK 

10. 

240.  0. 

000 

240. 

0. 

ODO 

15. 

150.  0. 

000 

240. 

0. 

000 

20. 

110.  0. 

002 

240. 

0. 

004 

25. 

75.  0. 

008 

240. 

0. 

<^5 

30. 

45.  0. 

010 

80. 

0. 

018 

35. 

25.  0. 

009 

25. 

0. 

009 

40. 

10.  0. 

006 

15. 

0. 

C09 

50. 

0.  0. 

006 

10. 

0- 

013 

25-MAY-S8 
15:  52:  40 
N  1  3  1. 
1 
7 

3 

4 

5 


NUMBER  OF  SYMPTOMS  -  44. 

SYMPTOM  CATEGORY  =  2 

NUMBER  OF  EXPOSURES  =  773. 

fraction  with  SYMPTOMS  =  0.05  7 
MODEL  INDEPENDENT  LIKELIHOOD  - 

TIME  AT  INDICATED  DEPTH 


-14.  301 3067 

at  fractional  risk  below  line 


depth  0.01000  0.02000  0.03000  0.04000  0.05000  0 


06000  0.  07000  0.  08000 


10.  00  -16.  1  -32.  3 

15.  00  -30.  6  -61.  5 

20.00  523.2  1051.8 

25.  00  24.  1  48.  5 

30.  00  10.  4  21.0 

35.  00  6.  1  12.  2 

40.  00  3.  7  7.  a 

45.  00  2.  6  5.  1 

50.  00  1.7  3.  4 

55.  00  1.1  2.  3 

60.  00  0.  8  1.6 

DOSE(F)  0.02216  0.04455 


-4S.  6  -65.  2  -81.  9 

-92.  a  -124.  3  -156.  2 

1585.7  2125.2  2670.4 

73.  2  78.  0  123.  2 

31.  6  42.  4  53.  3 

18.  4  24.  6  30.  7 

11.  7  15.  7  17.  7 

7.  7  10.  4  13.  0 

5.  1  6.  8  8.  6 

3.  4  4.  6  5.  7 

2.  4  3.  2  4.  0 

0.067  16  0.07001  0.11310 


-78.8  -115.7  -133.2 

-188.5  -221.0  -254.0 

3221.3  3778.1  4341.0 

148.  6  174.  3  200.  3 

64.  3  75.  4  86.  6 

37.  3  43.  8  50.  3 

23.  8  27.  9  32.  1 

1 5.  7  18-  4  21.2 

10.  4  12.  1  14.  0 

6.  9  8.  1  7.  3 

4.  8  5.  6  6.  5 

0.  13643  0.  16002  0.  18385 


DEPTH 
10. 
15. 
20. 
25.  ■ 
30. 
35. 
40. 
50. 


U.  S. 


OLD 

LIMI 

:ts 

TIME 

RISK 

240. 

0. 

000 

150. 

0. 

000 

110. 

0. 

002 

75. 

0. 

031 

45. 

0. 

042 

25. 

0. 

041 

10. 

0. 

026 

0. 

0. 

026 

NAVY 


NEW 

LIMITS 

TIME 

RISK 

240. 

0.  000 

240. 

0.  COO 

240. 

0.  005 

240. 

0.  095 

80. 

0.  074 

25. 

0.  041 

15. 

0.  038 

10. 

0.  058 

APPENDIX  C  -  SYMPTOM  CATEGORY  THREE 


,  25--MAY-88 

15:  34:  51 
N  1  3  3. 

1  1.  7313301563  = 

->  16.  1001243591 

3  0.  0840160698  =  Cl 

4  0.  1119670719  =  C2/C1 

Q  -291  48  =  MAXIMUM  LOG  LIKELIHOOD 


NUMBER  OF  SYMPTOMS  =  81. 

SYMPTOM  CATEGORY  =  3 

NUMBER  OF  EXPOSURES  =  773. 

FRACTION  WITH  SYMPTOMS  =  0. 10 5 

MODEL  INDEPENDENT  LIKELIHOOD  =  -24.9281578 


TIME  AT  INDICATED  DEPTH  FRACTIONAL  RISK  BELOW  LINE 
DEPTH  0.01000  0.02000  0.0300  0  0.04000  0.05000  0.06000  0.07000  0.08000 


10. 

00 

-23. 

3 

-46. 

.  9 

15. 

00 

-151. 

3 

-304. 

.  1 

20. 

00 

28. 

0 

56 

.  4 

25. 

00 

11. 

9 

23. 

.  9 

30. 

00 

6. 

7 

13 

.  4 

35. 

00 

4. 

3 

8 

.  6 

40. 

00 

2. 

9 

5 

.  7 

45. 

00 

2. 

0 

3 

.  9 

50. 

00 

1. 

3 

2 

.  6 

55. 

00 

0. 

9 

1 

.  8 

60. 

00 

0. 

6 

1 

.  3 

DOSE  (F) 

0.  01831 

0.  03^ 

68 

-70.  8 

-94.  8 

-119.  1 

453.  4 

-614.  4 

-772.  0 

85.  0 

113.  9 

143.  2 

36.  1 

48.  4 

60.  8 

20.  2 

27.  0 

33.  9 

12.  9 

17.  3 

21.  7 

a.  7 

11.  6 

14.  6 

5.  9 

7.  9 

10.  0 

4.  0 

5.  3 

6.  7 

2.  7 

3.  6 

4.  5 

1.  9 

2.  5 

3.  2 

055  49  0.  07437  0.  09345 


143.  7 

-168.  6 

-193.  7 

931.  3 

-1092.  2 

-1254.  9 

172.  7 

202.  5 

232.  7 

73.  3 

86.  0 

98.  8 

40.  9 

48.  0 

55.  2 

26.  2 

30.  7 

35.  3 

17.  6 

20.  7 

23.  7 

12.  0 

14.  1 

16.  2 

8.  1 

9.  5 

10.  9 

5.  5 

6.  4 

7.  4 

3.  a 

4.  5 

5.  2 

.  11273  0.  13221  0.  15191 


U.  S.  NAVY 

OLD  LIMITS  NEW  LIMITS 

DEPTH  TIME  RISK  TIME  RISK 

10.  240.  O.  000  240.  O.  COO 

15.  150.  0.000  240.  0.000 

20.  110.  0.039  240.  0.082 

25.  75.  0.  061  240.  0.  183 

30.  45.  0.  066  80.  0.  114 

35.  25.  0.  057  25.  0.  057 

40.  10.  0.035  15.  0.051 

50.  0.  0.  035  10.  0.  074 


